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The spectrum of meson and diquark excitations of dense color neutral cold quark matter is 
investigated in the framework of a 2-flavored Nambu-Jona-Lasinio type model, including a quark 
H- and color fj,s chemical potential. It was found out that in the color superconducting (2SC) 
. phase, i.e. at > /Xc = 342 MeV, /.ig aquires rather small values ~ 10 MeV in order to ensure 

the color neutrality. In this phase the tt- and a meson masses are evaluated around ~ 330 MeV. 
The spectrum of scalar diquarks in the color neutral 2SC phase consists of a heavy (SUc(2)-singlet) 
resonance with mass ~ 1100 MeV, four light diquarks with mass 3|/i8|, and one Nambu -Goldstone 
^ l' boson which is in accordance with the Goldstone theorem. Moreover, in the 2SC phase there are 

' five light stable particles as well as a heavy resonance in the spectrum of pseudo-scalar diquarks. 

, In the color symmetric phase, i.e. for /x < /ic, a mass splitting of scalar diquarks and antidiquarks 

is shown to arise it n ^ 0, contrary to the case of = 0, where the masses of scalar antidiquarks 
and diquarks are degenerate at the value ~ 700 MeV. If the coupling strength in the pseudo-scalar 
diquark channel is the same as in the scalar diquark one (as for QCD-inspired NJL models), then 
^>f-^ ■ in the color symmetric phase pseudo-scalar diquarks are not allowed to exist as stable particles. 

> . 
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■ I. INTRODUCTION 

Dh: . . . 

I Recent investigations, performed in the framework of perturbative QCD, show that at low temperatures and asymp- 

Oh! toticaUy high values of the quark chemical potential /i the dense baryonic matter is a color superconductor Ev- 
^ idently, at rather small values of /z a more adequate description of this phenomenon can be done with the help of 
• • . different effective models, such as Nambu - Jona-Lasinio- (NJL) type field theories with four-fermionic interaction 
' 0i -'-'^ such a way, on the basis of NJL- type models with two quark flavors it was shown that the color super- 
k>( I conducting (2SC) phase might be yet present at rather small values of ~ 350 MeV, i.e. at baryon densities only 

■ several times larger, than the density of ordinary nuclear matter (see reviews |1,Q). (For simplicity, throughout the 
paper we will study quark matter, composed from two quark flavors, i.e. up- and down quarks, only.) This is just the 
density of compact star cores. So, color superconductivity, which is eventually existing inside compact stars, might 
influence different observable astrophysical processes, and thus deserves to be studied in more details. 

In the early NJL approach to color superconductivity |^ 0| , the density of the color charge Qg was not constraint 
to zero in the 2SC ground state (the densities of other color charges Qi, where i = 1, .., 7, are zero in this phase J\) 
leading to a nonvanishing difference between the densities of red/green quarks and blue ones. In this case, the mesonic 
and diquark excitations of dense quark matter were considered in the framework of a simple two-flavored NJL model 
with a single quark chemical potential In particular, it was shown that in the color asymmetric 2SC phase 

of this model does arise an abnormal number of three Nambu - Goldstone bosons (NG) instead of expected five ^ , 
and TT-mesons are stable excitations of its ground state with masses ^ 300 MeV. Besides, there are two light stable 
scalar diquark modes, whose masses are proportional to (Qs), as well as one heavy scalar diquark resonance in the 
2SC phase. 

In reality, however, i.e. possibly in compact star cores or in relativistic heavy ion experiments, there are several 
physical constraints on the quark matter. The most evident one is its color neutrality, which means a vanishing of 
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^ Recall, that in the 2SC phase of the 2-flavored NJL model the initial SUc(3) color symmetry is spontaneously broken down to the 
SUc(2) one, so one might expect naively five massless bosons. However, due to Lorentz noninvariance of the system, there are indeed 
only three NG bosons. 
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a bulk Qi color charges {i = 1, ...,8). Indeed, since the lump of quark matter, which might be created after heavy 
ion collisions, originated from color neutral objects and is surrounded by a color neutral medium, it is expected 
to be globally color neutral. To fulfill this requirement, usually the local color neutrality constraint is imposed by 
introducing several new chemical potentials, /i3, /ig, etc, into a NJL model 0, 0, ^|. Otherwise, there will be 
produced a chromo-electric field resulting in the flow of color charges, so a homogeneous and color conducting quark 
medium with nonzero color charge densities is not allowed to be a stable one (Recall, there is no need to add 

new chemical potentials into QCD. The point is that in the QCD 2SC ground state a nonzero value of the eightth 
gluon field component might be generated. Effectively, it is the /xg-chemical potential, so color neutrality is fulfilled 
automatically in the QCD approach [ij.) 

In the present paper, in contrast to our previous investigations 0,01' study the mesonic and diquark excitations 
of color neutral quark matter ^ in the framework of a simple two flavor NJL model at zero temperature. We consider 
both the case of rather small values of the quark chemical potential fi (the SUc{3) color symmetric (normal) phase), 
and the case of /x-values, corresponding to the 2SC phase of the model. In addition, the properties of pseudo-scalar 
diquarks are also included into the consideration. 

The paper is organized as follows. In Section II the thermodynamic potential as well as the effective action of 
the NJL model, extended by an additional color /^g-chemical potential term, is obtained in the one quark loop 
approximation. Further, in Section III, the gap equations and the phase structure of the model is investigated under 
the color neutrality constraint. Here the values of /ig are obtained at which the 2SC phase is a color neutral one. Then, 
in the Sections IV, V, and VI, the masses of the tt- and a mesons, scalar diquarks, and pseudo-scalar diquarks are 
considered, respectively. Finally, in the Appendix, the influence of the mixing between cr-meson and scalar diquarks 
on the (T-mass is discussed. 



II. THE MODEL AND ITS EFFECTIVE ACTION 

In the original version of the NJL model the four-fermionic interaction of a proton p and neutron n doublet 
was considered, and the principle of dynamical chiral symmetry breaking was demonstrated. Later, the {p, n)-doublet 
was replaced by a doublet of colored up u and down d quarks, in order to describe phenomenologically the physics of 
light mesons 

El III 111 111 , diquarks [H , as well as meson-baryon interactions 0, 12^ . In this sense the NJL 
model may be thought of as an effective theory for low energy QCD. (Of course, it is necessary to remember that 
in the NJL model, in contrast with QCD, quarks are not confined in the hardronic phase, which is a shortcoming 
of the model.) At present time the phenomenon of dynamical (chiral) symmetry breaking is one of the cornerstones 
of modern physics. So, this effect was studied in the framework of NJL-type models under the influence of external 
magnetic fields in curved space-times |^, in spaces with nontrivial topology jH, etc. Formally, as it was 
mentioned above, quarks are presented in the mass spectrum of the model. So it is very suitable for the description 
of normal hot and/or dense quark matter |l5l l24l I25ll26| | in which, as it is known, quarks are not confined. NJL-type 
models still remain a simple but useful instrument for studying color superconducting quark matter at moderate 
densities H^IEHi where analytical and lattice computations in the framework of QCD are impossible. 

We start from the following 2-flavor NJL Lagrangian, called for the description of interactions in the quark-antiquark, 
scalar diquark- as well as pseudo-scalar diquark channels at low and moderate energies and baryon densities (the 
consideration is performed in Minkovski space-time notation): 

Lg ^ q^y^icl^, ~mo + fiJ° q + G\^{qqf + {qij'^fqf} + 
+ {Hs[fil'^T2XAqWl'^T2\Aq^] + Hp[q'^T2\Aq][qT2\Aq'^]}, (1) 

A=2,5,7 

where > is the quark chemical potential, the quark field q = qia is a flavor doublet and color triplet as well as 
a four-component Dirac spinor, where z = 1,2 (or i = u,d) and a = 1,2,3 (or a = r,g,b). q^ = Cg*, (f = q*C 
are charge-conjugated spinors, and C — 47^7° is the charge conjugation matrix (symbol t denotes the transposition 
operation). It is supposed that up and down quarks have equal current (bare) mass toq. Furthermore, we use the 
notations Ta for Pauli matrices, and Xa for antisymmetric Gell-Mann matrices in flavor and color space, respectively. 



^ Note that in compact star cores one should consider the electrically neutral quark matter in beta equilibrium, whereas in heavy-ion 
experiments the isospin and strangeness are conserved quantities. For simplicity , these additional constraints on the quark matter are 
ignored in the present consideration. 

^ Indeed, consider two-flavor QCD, symmetric under the color group SUc{3). By integrating in the generating functional of QCD over 
gluons and further "approximating" the nonperturbative gluon propagator by a 5— function, one arrives at an effective local chiral 
four-quark interaction of the NJL type describing low-energy hadron physics. Moreover, by performing a Fierz transformation of the 
interaction terms, it possible to obtain a NJL-type Lagrangian describing the interaction of quarks in the scalar and pseudo-scalar {qq)- 
as well as scalar and pseudo-scalar diquark (qq) channels (see, e.g., the Lagrangian below). 
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Clearly, the Lagrangian Lq is invariant under transformations from the color SUc(3)- as well as baryon Ub(1) groups. 
In addition, at mg = this Lagrangian is symmetric under the chiral SU(2)l x SU(2)r group (chiral transformations 
act on the flavor indices of quark fields only). Moreover, since Q — + B /2, where = T3/2 is the generator of 
the third isospin component, Q is the generator of the electric charge, and B is the baryon charge generator, our 
system is symmetric under the electromagnetic group Uq(1) as well. If the Lagrangian |^ is obtained from the QCD 
one-gluon exchange approximation, then G : Hg : Hp = 4:3:3. However, in the present consideration we will not 
fix relations between coupling constants, so they are treated as free parameters. It is necessary to note also that at 
/X = the Lagrangian is invariant under the charge conjugation symmetry (g — > q'^ = Cg*, q ^ (f = q*C) that 
is, however, spoiled by the chemical potential term. 

Furtheron, the temperature is chosen to be zero throughout the paper. In this case there is a critical value /ic of the 
chemical potential, such that ai fi < jlc the color symmetric (normal) phase is realized in the system (evidently, the 
ground state of this phase is a color singlet). However, at fj, > flc one obtains a color superconducting (2SC) phase in 
which SUc(3) is spontaneously broken down to SUc(2). Only two quark colors, say red and green, participate in the 
gap formation in the 2SC phase, the blue quarks stay gapless. So the densities of red/green quarks, Ur^g, are equal 
in this phase, however the density Ub of blue quarks is not equal to n,-.g, i.e. local color neutrality is broken. ^ To 
restore local color neutrality in the 2SC phase of the model ([Q, usually an additional chemical potential term fisQs is 
introduced into the considerations Q, where Qs = q^^T^q, and Tg =diag(l, 1, —2) = %/3A8 is the matrix in the color 
space. Hence, 

Lg ^ L = Lq + ^sQs- (2) 

If /ig is an independent model parameter, then Lagrangian |(2Jl is symmetric under the color SUc(2) x Ua8(1) group. 
However, if local color neutrality is imposed, then the chemical potential fig, is not an independent parameter. Its 
value must be chosen in such a way that the ground state expectation value (Qs) is identically equal to zero. Hence, 
fis depends on /i etc. For example, in the color symmetric phase, i.e. at < /ic (it will be shown in section UTTl that 
in general case /ic flc), where (Qs) = even in the theory we have to put /ig = 0. However /ig has a nontrivial 
/i-dependence at fj, > fic in order to supply the zero value of (Qs)- As a consequence, we see that the color symmetry 
group of the model Q depends on fi: at /i < fic it is SUc(3), whereas at ^ > it is the SUc(2) x Ua8(1) subgroup 
of SUc(3) which is just the color symmetry group of the term /zsQs- 

In the present paper we are going to study both the ground state properties of the system with Lagrangian L, and 
its mass spectrum in the quark-, meson- and diquark sectors. So, we have to obtain the thermodynamic potential 
as well as the effective action of the model up to a second order in the bosonic degrees of freedom. To begin with, let 
us introduce the linearized version of Lagrangian L that contains auxiliary bosonic fields: 



L 



Yid„ + fl^ + l-lsTg - (T - Too - 47 TTaTa 



AS* AS A P* A, P 



4Hs ^ ^ 4Hp 



- ^[ifil'r^XAq] - ^[q^J'T,XAq''] - ^[q^r^XA'q] ~ ^[qr,XA'q% (3) 

where as well as in the following the summation over repeated indices a = 1,2,3 and A, A' = 2,5,7 is implied. 
Lagrangians L and L are equivalent on the equations of motion for bosonic fields, from which it follows that 

aix) ^~2Giqq), /\\{x)^-2Hs{fl-i''T2XAq). A'X{x)=^-2Hsiqij''T2XAq^), 
TTa{x)^-2G{qi-i\aq), A\,{x) ^ -2Hp{f t^Xa' q) , AP;,ix) ^ ~2HpiqT2XA'q^). (4) 

One can easily see from Q that mesonic fields cr, tTq are real quantities, i.e. {a{x)y — (j{x), {7Ta{x)y — Tra{x) (the 
superscript symbol f denotes the hcrmitian conjugation), but all diquark fields A^^ are complex ones, so 

(A^(x))t = Ai*(x), (AP,(x))t = AP:(x). 

Moreover, A^ and A^, are scalars and pseudo-scalars, correspondingly. Clearly, the real cr and tTq, fields are color 
singlets, all scalar diquarks A^ form a color antitriplet 3c of the SUc(3) group. The same is true for pseudo-scalar 
diquarks A^, which are also the components of an 3c-multiplet of the color group. Evidently, in the SUc(3)-color 
symmetric phase (at fj, < /Zc) all diquark fields must have zero ground state expectation values, i. e. (A^) — and 
(A^,) = 0. Otherwise, we have an indication that the ground state of the system is no more an SUc(3)-invariant one. 

Lagrangian Q provides us with a common footing for obtaining both the thermodynamic potential (TDP) and the 
mass spectrum for bosonic excitations. Indeed, in the one-fermion loop approximation, the effective action for the 



* Thus, (Qs) 7^ 0- However, the color charge Qg = q'y^T^q, where T3 =diag(l, —1,0) is the matrix in the color space, vanishes in the 
2SC phase. Other color charges Qi [i 7^ 8) are also zero in this phase 0. 
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boson fields is expressed through the path integral over quark fields: 

exp(z5eff(a,7ra, A^^, A^f*)) = N' [dq][dq]e^p(^i 



Ld^x 



where 



d^x 



AG AHs AHp 

N' is a normalization constant. The quark contribution term Scs is here given by: 



ScS, 



exp(i5eff) ^ N' [dq][dq]exp(^^ 



qD+q + fD-q^ - qKq^ - f K* q 



d^x 



(5) 



(6) 



In © we have used the following notations 

£)+ = ij^d,, - mo + - S, D~ = i^^d^ - mo - fij° - S*, S = cr + ij^naTa, 
E* = a + i^\arl K* = (A^* + zAX7')t2Aa, ^ = (A^ + *A^75)r2A^, 



(7) 



where are nontrivial operators in the coordinate, spinor, color and flavor spaces.'"^ In the framework of the Nambu- 
Gorkov formalism, where quarks are composed into a bispinor ^ = ^ ^ , it is possible to integrate in ^ over quark 



fields and obtain 



5off(cr,7rQ, A^^, A^D = ^Tr{jvGs/cx} In 



D+ , -K 
-K* D- 



(8) 



Besides of an evident trace in the two-dimensional Nambu-Gorkov space (NG), the Tr-operation in (jHJ stands for 
calculating the trace in spinor- (s), flavor- (/), color- (c) as well as four-dimensional coordinate- (x) spaces, corre- 
spondingly. 

Starting from (0)-®, it is possible to define the thermodynamic potential ri((T, tt^, A^^, A^f*) of the model: 



5, 



•off 



CT.TTa ,A^.'^ ,A^'f* —const 



(9) 



where, in the spirit of the mean-field approximation, all boson fields are supposed to be x-independent. It is well- 
known that ground state expectation values {cr{x)) = ct°, {TTa{x)) = 7r°, (A^^(x)) = A^^°, (A^f*(a;)) = A^f*° are 
coordinates of the global minimum point of the thermodynamic potential f2, i.e. they form a solution of the gap 
equations (evidently, in our approach all ground state expectation values do not depend on coordinates x): 



0, 



dfl 
dAf* 



0. 



(10) 



Let us make the following shift of bosonic fields in iSeff: cr(x) a{x) + a°, Tra{x) T^a{x) + 7r°, A^^*(x) 
A^^*(x) -I- A^^*°, A^^(x) A^^(x) -I- A^^°, where a", 7r°, A^^°, A^^*° have no coordinate dependency. In this case 
the matrix Z from (|HJ| is transformed in the following way: 



-k; , D- 



E , K 



- c-l 




K* , S 

where 5*0 is the quark propagator matrix in the Nambu-Gorkov representation, and 

{Ko, K:,Dt, Eo, E*) = (K, K*,D^, E, E* 



E , K 



(11) 



(T—a° .7r„ — 7r° 



Then, expanding the expression JSJ up to a second order over the meson and diquark fields, we have 

5eff (a, TT,, A^, A^D = Sj^'j + Sj^'j (a, ^„ A^^ A^f*) + • • • , 



(12) 



^ In order to bring the quark seetor of the Lagrangian to the expression, given in the square brackets of JHJ, we have also used the 
following well-known relations: 9' = -d^, CYC'^ = -(YY, C-y^C''^ = {-f^Y = 7^, r^fr^ = -(f)*, = ( ° ' ~M . 
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where (due to the gap equations, the term hnear over meson and diquark fields is absent in (|12|l ') 



c(0) _ 
'-'eff - ^ 



4G 



A"-^A ' 

4i?„ 



;Trj7VGsc/a;} In [Sq ^) — 



(13) 



5,^^^(a,^,,A7,Ar) = - 



^2+^2 ^P^^P* 



4G 



4iJ., 



4ir„ 



:(2) 



(14) 



In the foUowing on the basis of the effective action we wih study the spectrum of meson/diquark excitations in 
the color superconducting phase of the NJL model under consideration. So, it is convenient to present the effective 
action H14|l in the following form: 



where 



5(2) = - 



c(2) ^ <-(2) , <^(2) (,(2) 

'-^off '-'mesons '-^diquarks '-^mixcd' 



S. 



(2) 

diquarks 



d^x 



A^Af 



-Tr,,/, {5i2if*5i2if* + 2S11KS22K* + S21KS21K} , 



c(2) ^ 1 
mixed 2 



4ii, 4iip 

Tr,,/, {5iil]5i2if* + S2i^SiiK + Si2^'S22K* + S2iKS22^'} 



(15) 

(16) 
(17) 
(18) 



and Sij are the matrix elements of the quark propagator matrix 5*0, given in ()ll|l . 



III. GAP EQUATIONS AND COLOR NEUTRALITY CONDITION 



Let us for a moment assume that in ((5J) the chemical potential fig is an independent parameter (^ 0). Then, 
Lagrangian L is invariant under the color SUc(2) x Ua8(1) symmetry group. Recall, that the phase structure of 
any theory is defined by a competition of its order parameters. In our case, the order parameters (ground state 
expectation values), (cr(a;)), {TTa{x)), (A^(a;)), (A^,(a:)), are obtained from a solution of the gap equations (see the 
previous section). Since the consideration of the model l(2Jl with a total set of order parameters is a very hard task, 
we shall assume that parity is conserved, i.e. {iraix)) = 0, (A^, (x)) = (in the section IVl Al some arguments are 
presented, however, that at sufficiently high values of pseudo-scalar coupling Hp parity might be spontaneously broken 
down), thus having to deal only with (cr(x)) and (A^(x)). In this case three different phases might exist in the model 
©I i) In the first one, the normal phase, (A^(x)) = for all A = 2,5,7. In this phase the initial color symmetry 
remains intact, ii) The second one is a well-known 2SC phase with (A2(a::)) ^ and (A^ j{x)) — 0. The ground state 
of this phase is invariant under SUc(2)-color symmetry, iii) Finally, there might exist a phase with (A|(a;)) ^ 0, 
(A|(a;)) ^ 0, and (Af (a:)) = 0. (Note, the two phases ii) and iii) are not unitarily equivalent, since there are no color 
transformations from SUc(2) x UAg(l) that connects the corresponding ground state expectation values.) However, 
since color neutrality cannot be achieved in the ground state of the form iii) (see [tI), throughout of our paper only 
two order parameters, {^{x)) and (A|(a;)) = A, will be taken into account whereas other ones will be supposed to have 
zero expectation values, i.e. {TTa{x)) — 0, (A^, (a;)) ~ 0, (Ag -^(x)) — 0. So, below only the competition between the 
normal phase ( A = 0) and the 2SC one ( A ^ 0) will be considered. As a result, one may deal with a thermodynamic 
potential Q which depends on two variables A, (a) or, equivalently, A,m = toq + (f), only. Then, the expression for 
the thermodynamic potential Q can be calculated with the help of (|13|l (sec also, e.g., [ij): 



n{ni,A) = 



(m — TOo)^ 
4G 



'iH, 



± 



_fq_ 



\E 



A I 



2E 



d^q 
(2^)3 



(19) 



where = ^J{E^y + jAp, £'± = £^±/I, E^ = E±il, p, = ^ + fj,^^ fl = ^J.~2fis, E = ^ ^ w? . Since the integrals 



in the right hand side of (|19|l are ultraviolet divergent, we regularize them and the other divergent integrals below by 
implementing a three-dimensional cutoff A. The resulting gap equations look like: 

4iA / . ^ I ' I r, A / y 



(2n) 



m — mo 
2G 



Aim J 



d*q 1 r E^ 
J2^E\J2 



1 „ ^ f d'^q I ( E^ 



(20) 



}■ 



(21) 
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In (|20() . H21() as well as in other expressions, containing four-dimensional momentum integrals, qo is shorthand for 
Qq + ie ■ sign(go), where e — *■ 0_|_. This prescription correctly implements the roles of /i, fig, as chemical potentials and 
preserves the causality of the theory (see, e.g. p^l. 

Now, let us impose the local color neutrality requirement on the ground state of the model. It means that the 
quantity /ig takes such values that the density of the 8-color charge (Qs) = —dQ/dfis equals zero for arbitrary fixed 
values of other model parameters. So, we have the following local color neutrality constraint: 

dn f (fq cE+ E- E+ E- ^ 

The system of eqs. H2()(l - (l22ll has two different solutions. As we have already discussed after (O, the first one (with 
A = 0, /i8 = 0) corresponds to the SU(3)c-symmetric phase of the model (normal phase), the second one (with A ^ 
and ^8 7^ 0) to the 2SC phase. As usual, solutions of these equations give local extrema of the thermodynamic 
potential ri(TO, A) (|19|) . so one should also check which of them corresponds to the absolute minimum of O. Having 
found the solution corresponding to the stable state of quark matter (the absolute minimum of we obtained the 
behavior of the gaps m, A as well as the fig vs. the quark chemical potential /i (see Figs 1,2). (Note, in all numerical 
calculations of the paper we use the parameter set 

G = 5.86GeV~^ A = 618 MeV, toq = 5.67 MeV, i?, = 3G/4 (23) 

that leads in the framework of the NJL model to the well-known vacuum phenomenological values of the pion weak- 
decay constant = 92.4 MeV, pion mass A/^r = 140 MeV, and chiral quark condensate {qq) = —(245 MeV)'^. The 
relation between Hs and G in (|23|l is induced, e.g., by the structure of QCD four-quark vertices in the one-gluon 
exchange approximation.) The region fi < fi^ = 342 MeV is the domain of SUc(3)-color symmetric quark matter 
because £7 in this case is minimized by m ^ and A = 0, fig = 0. For fi > fi^, the solution with m 7^ 0, A 7^ and 
fig ^ 0, corresponding to the 2SC phase, gives the global minimum of fi, and thereby the color superconducting phase 
is favored. The transition between these two phases is of the first-order, which is characterized by a discontinuity in 
the behavior of m and A vs. ^ (see Fig. 1). Finally, remark that the critical value /Etc of the quark chemical potential 
was calculated in the framework of the model as well. In terms of the parameter set (I23|l we have jlc = 350 MeV 
0,0], i.e. in the color neutral quark matter the transition to the 2SC phase occurred at slightly less values of the 
quark chemical potentials. 

Similarly to (2Si] , it is possible to find the following expressions for the matrix elements Sij of the quark propagator 
matrix Sq: 

A. /^%_^^|^(.)^ (24) 



(27r)4 [qa+E+ qo ~ E' 



"~ '(2.)^^ U-{E-^r^''^^qi-iEir^''-i'^'' + 

d'l (^:!A^ + -^}pt\ (25) 



(27r)4 [qn+E- qo - E 



S,, = -^A*r,X, / 7:^e--(-^) + + \ , (26) 



S.. - -A^.A. / ^e-^(^-) + -,2^ } , (27) 



(2^)4 [q^-iE^r ql'iElY 

where A± = i(l ± 2Lil|z!!i))^ and P^^ — diag(l, 1,0), P^'^^ ~ diag(0,0, 1) are the projectors on the red/green and 
blue directions in the color space, correspondingly. 

The poles of the matrix elements (|24|) - (|27f) of the quark propagator give the dispersion laws, i.e. the momentum 
dependence of energy, for quarks in a medium. Thus, we have E^ for the energy of red/green quarks and E'^ for the 
energy of red/green antiquarks. Moreover, the energy of blue quarks (antiquarks) is equal to E^ (E^)- It is clear 
from Figs 1,2 that in the 2SC phase, i.e. at /i > /Zc, we have {fi + fig) > m, {fi — 2fig) > m. Then, in this phase the 
quantity E can reach both the value {fi + fig) (the Fermi energy for red/green quarks) and the value {fi — 2fig) (the 
Fermi energy for blue quarks). In this case, in order to create a red/green quark in the 2SC phase, a minimal amount 
of energy (the gap) equal to |A| at the Fermi level {E = fi + fig) is required. Similarly, there is no energy cost to 
create a blue quark at its Fermi level E = fi — 2fig, i.e. blue quarks are gapless in the 2SC phase. In the normal phase 
of the model, i.e. at fi < fic, the minimal energy (to — fi) is needed for the creation of a quark of any color. Note, 
both in the 2SC and normal phases of the model the minimal energy required for a quark creation differs from the 
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FIG. 1: The constituent quark mass m vs /i and FIG. 2: The behaviour of fis vs /i in the model 

color gap A vs fi under local color neutrality con- JJJ when local color neutrality is imposed, 

straint. 



minimal energy required for the creation of an antiquark of the same color. This fact reflects the breaking of charge 
conjugation symmetry in the presence of chemical potentials. 

Without loss of generality, we will assume throughout the paper that A is a real nonnegative quantity. Given 
the explicit expression for the quark propagator S'o, in the next sections we will calculate two-point (unnormalized) 
correlators of meson and diquark fluctuations over the ground state in the one-loop (mean-field) approximation and 
find their masses. 



IV. MASSES OF THE tt- AND a-MESONS 

It turns out that the S^^^^^^ -part ifT^ of the effective action is composed from a{x), A|(a;) and A|*(a;) fields only. 
So it provides us with nondiagonal matrix elements T^x = A|*, A|) of the inverse propagator matrix of cr, A|*, 

(2) 

and Aj. Moreover, each term in iS^/xcd proportional to A or A* as well as to the constituent quark mass m (see 
Appendix). Hence, in the color symmetric phase (A = 0) there is no mixing between cr-meson and diquark A| at 
all. The parameter m is small (or even equals zero if mo — 0) in the 2SC phase, so we ignore for simplicity the (7-A2 
mixing effect in this phase, too. As a result, in order to get the masses of mesons we use only the effective action (|16f) 

which has the form 5mcsons = S^V + S^J , where 

^2 



/G i 
+ -Tr^c/a; {SiidSiid + 2Si2aS2icr + S22<jS22<j} 



= -i y" d^ud^v (j{u)T{u - w)cr(w), (28) 

S^^} = - y d^X^ + '-T^scSx {5ll(*7'^aTa)5ii(l7'7r,Tb) + 2^12 (it'tT.T* ) ^21 (*7'Tf,rfa) + 

+ 522(«7'^ar*)522(*7'7rbT*)} = j d^ud\ TTk{u)Iiu{u ~ v)ni{v) . (29) 

In these formulae T{x — y) is the inverse propagator of (T-mesons, and ^ab{x — y) is the (diagonal) matrix of the inverse 
TT-meson propagator. Evidently, that 

J2 ^2 ^^2) 

T{x-v)^- nab{x-y) = - I I \ - 

da[y)b(T{x) dnbiyjOTTaix) 

Next, using in H28|l - H29l) the expressions (|24|l - 127|l for the matrix elements Sij, it is possible to obtain with the help 
of relations the functions r(a:; — y), Ilab{x — y), and then their momentum space representations, r(p), I{ab{p), 
correspondingly.^ The zeros of these functions determine the particle and antiparticle dispersion laws, i.e. the relations 



One should not be confused by the coordinate or momentum space representation which is used for the inverse propagators and other 
quantities, since this is clear from the arguments of these functions or the context. 
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between their energy and three-momenta. In the present paper, we are mainly interested in the investigation of the 
modification of meson and diquark masses in dense and cold color neutral quark matter. In this case, the particle 
mass is defined as the value of its energy in the rest frame, p = (see, e.g., |3,|23)i where the calculation of inverse 
propagators for a- and 7r-mesons is significantly simplified. Indeed, in the rest frame it is possible to get: 



2G ""7 (2^)3 E+E^ ^e++E^)^-pI 



(Pq d{E~il)E 
(2^ AE^-pI 



-16<5a6 / ,._v. 2 = SabTliPo), (31) 



T{po)=To{pl)+T,{pl), (32) 
1 ^ f d^q f E^E^ + E+ E- + A"^ E^ + E^ f (fq 9{E - fl) 



ro(Po) - 2G V (2^)3 E^ E+E^ {E+ + E^^ - pi J i2n f E AE^ pi ' ^^^^ 



J {2^YE^\Ei[pl-A{EX?] E-^[pl-A{E-^Y]^' ^' 

where the same notations as in (|19f) were used. The zeros of the functions (|31|l . (|32|l give us the masses of tt- and 
cr-mesons, respectively (see Fig. 3). ^ We see that in the 2SC phase the masses of the sigma- and pi-meson are about 
300 MeV. Moreover, the pion is a stable particle in this phase (only electroweak decay channels are allowed). This 
conclusion is supported by the following arguments. It is clear that the first and the second integrals in (|31l) are 
analytical functions in the whole complex Pq plane, except the cuts Ef^-^^^ < pi < oo and (2/1)^ < pi < oo, respectively 
(here Emm = y/ {ft — 'niY + | Ap -I- \/{ft + m)^ + | is the minimum of the expression E^ + E'^, which is taken at 
the point \p\ — 0)- Evidently, Emin corresponds to the threshold for the pion decay into a red-green quark- antiquark 
pair, whereas 2/i corresponds to the threshold for the pion decay into a blue quark-antiquark pair. It is easily seen 
from Fig. 3 that in the 2SC phase the pion mass is less than the values of these two thresholds. Because there are 
no other singularities in H31|l corresponding to different channels of the pion decay, we can conclude that in the 2SC 
phase the pion is a stable particle. 

The neglection of the mixing between the a-meson and A2-scalar diquark also results in a stable tr-meson. However, 
if the mixing is taken into account, then in the 2SC phase the cr-meson is a resonance, decaying into a pair of quarks, 
whose width is a rather small quantity, i.e. about 30 MeV (see Appendix). 

V. MASSES OF SCALAR DIQUARKS 

As in the previous section, we will ignor for simplicity the term S^^j^^^ (|18|) . which mixes a and A|*, Aj diquarks, in 
the effective action (|15|) . In this case, in order to obtain the masses of diquarks, we need to analyze the term 5^^^^^^^^.]^^ 
()17|l only. It can be easily presented in the following form 

'^diquarks X! { '^i ^1 A + '^p^ A A } ' (35) 

A=2,5,7 

where labels s,p denote the contributions from scalar- and pseudo-scalar diquark fields, correspondingly, and 

S^\j, = -J d^x^^^ + ^Tr,,/, {^n*A^7V2AA522«AX7'T2A^} , (36) 
4'aa = - / d'x^f^ + '-Tr^^f. {^n A^^r2A^^22 APV2A^} , (37) 

for fixed A = 5,7, and 

= -J d^^^^ + '-Trscf. {^i2zAr7'^2A25i2iAf 7V2A2+ 



In our numerical investigations of the 2SC phase, presented in Fig. 3, we have ignored in I32i the term ri(pg) proportional to A^m,'^, 
since it is comparable (or even less) in magnitude with nondiagonal elements Ta-xipo) = -^2* '^2) '^^ ^^'^ ^'^^^ inverse propagator 
and Aj (see Appendix), which are not taken into account in the above consideration. 
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FIG. 3: The masses of the cr-meson and pion as functions of /i, when mixing of a and A2 is neglected. 



2^nzA^7V2A2^22zAr7^T-2A2 + ^2liA27'^T-2A2^21»A27^T2A2} , 



(38) 



^(2) 
'p 22 



PAP* 



4Hp ' 4 



TTr^c/a; {5i2A^*r2A2S'l2Af*T2A2 



+ 25nA^T2A2522Af T2A2 + 52iAfT2A252iA|r2A2} . (39) 

It follows from (|35|) - (|39|l that there is no mixing between scalar- and pseudo-scalar diquarks. Moreover, scalar diquarks 
(or pseudo-scalar ones), as such, are not mixed to one another. Starting from the above formulae, we will find the 
inverse propagators of diquarks which migth be introduced by the following way (A — 2,5,7): 



r(2) 



AA 



(2) 



p AA 



X.Y 



2 ^ 



d'^ud'^v X{u)T^%.{u - v)Y{v), 
d'^ud'^v P{u)T'^^q{u - v)Q{v), 



(40) 
(41) 



where (for each fixed value of yl) X{x),Y{x) = A%{x) , A"^ (x) , P{x),Q{x) = (x), A^* (x) and r^V(^) or ^pq{z) 
are matrix elements of the 2x2 inverse propagator matrix for A^(x), A^*(x)- or A^(a;), A^*(a;)-fLelds, respectively. 
Given diquark propagators, it is possible then to obtain the masses of diquarks. 



A. Scalar diquarks in the 2SC phase (A 7^ 0, /xg 7^ 0) 

In the present section, using (|40|l for different A — 2,b, 7, we will study step-by-step the masses of scalar diquark 
excitations in the 2SC phase of the model (O, when the color neutrality condition is taken into account. Let us begin 
with the A| -A|* diquark sector. It follows from at A = 5 that 



'5Y{y)5X{x) 



(42) 



(recall, X,Y = A|(a;),A|*(x) in the case under consideration). Note also that r^*(z) is a symmetric matrix, i.e 



p5s 
^ XY 



p5s / 



-z). It is clear from (gOl) and (g^l that T^^s^.{z) 



0, and this matrix has nonzero 



elements of the form (in the momentum space representation) : 



(43) 
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r%A'^>{p) = ^ -^Tr.^l ^ {522(g+p)zf A5^ii(g)zf As} , (44) 

where the Fourier-transformed expressions Sii{q), 5*22(9) can be easily determined from (|24() . (|25|l . It follows from 
(I43|l - H44|l that ^^Af'Af{~P) = r^s^j.(p). Since we are interested in diquark masses, it is necessary to use the rest 
frame in i.e. p — (po, 0,0,0) (see also 0,11^). In this case the calculation of matrix elements is 

greatly simplified, and mass excitations are connected with the zeros of the quantity det r^''(po) in the pp-plane.^ So, 
we have 



T^s , ^ 1 ^. / A qo + E+ qo-E- 

ArAim) V (27r)M(p„+^„+^+)(g2_(^+)2) (p„+^„_^-)(^2_(^-)2)J- 



The Po, go-dependency in the integrand of H45() is presented in an evident form. Other quantities in H45() . such as E^ 
etc, depend on |<7| only. The expression H45|) is valid both for A = and A ^ 0. For the case A ^ 0, i.e. in the color 
superconducting phase, it is possible to use the gap equation (|20|l in order to eliminate the coupling constant from 
this formula. In this way we find: 



rrrA;(Po) = 4z(po-3/i8) 



dS 1 



(27r)4 1^ + ^„ + E+)iq^, - (E+r) (po + go - E-){qi - (E^)^) j 
= 2(po - 3M8)ff(po), ri^iAr (Po) = -2(po + 3^M)H{-po). (46) 

The last equation in (|46|l is due to the relation ^^s*^f{—po) — rA=A"(Po)- Recall, in 14t)|l go + Po and go are 
shorthands for (po + qo) + • sign(po + 90) and qo + ie ■ sign(go), where e ^ 0+ (see also comment after formula (|21f) '). 
Performing in H46|l the go-integration, one obtains 



H{po) 



d^q j 9{E+) 0{-E+) 9{E- 



(2^)' i (po + E+ + E+)E+ (po +E+~ E+)E+ (po - E- - E^)E^ 

{po -E-+ E^)E^ J 



(47) 



Now, with the help of (|46|l one easily gets the expression for the determinant of the inverse propagator matrix r'^''(po) 
in the rest frame: 

detr5^(po) = T%.^.{po)T%^,^,{po) ^ -A{pI - 9fil)H{po)H{~po). (48) 

The diquark squared mass spectrum in the A|-sector of the theory is defined by zeros of the detr'^'*(po) in the p^- 
plane. Evidently, the point p^ — 9/z| is the solution of the equation detr^''(po) = 0. Let us suppose that some nonzero 
point Po — —Mfy is the zero of the function H{po), i.e. H{—Mfy) = 0. Then, at po — Mfy the determinant H48|l is 
also equal to zero, so the point Pq = (M|))^ is another zero of detr^^(po) in the pp-plane, and the second bosonic 
excitation of this sector has nonzero mass Mfj. It follows from H47|l and H22(l that H{po) is proportional to (Q8) at 
the point po = Sf^s- Namely, (S/xg) = ~{Qs) / (4A-^). Since (Qg) is zero due to the constraint we may conclude 
that Mfy = Sj^sl- Hence, in the A|-sector of the model there are two real bosonic excitations with equal masses 
MS = 3|/i8|. 

The similar is true for the A^-sector of the model, so in the whole A|, Ay-sector of the NJL model that is under the 
color neutrality constraint there are four massive scalar excitations with equal masses M|) = 3|//8|- These particles 
form two real antidoublets of the SUc(2)-group, or one complex antidoublct. 

Consider now the diquark excitations of the 2SC ground state in the A|-A|* sector of the model. In this case, the 
matrix T'^''{po) (the momentum representation for the inverse propagator matrix at p = 0, i.e. in the rest frame) has 
the following structure: 

ri'^iAi (po) = r'^r^r ^^o) = 4A2/o(p^), ri''.^„ (po) = r^^..^. (-po) = 

= (4A2-2p2)/o(p2) + 4po/i(p2), (49) 



Recently, the Bethe-Salpeter equation approach ha s be en used to obtain diquark masses in the 2SC phase of cold dense QCD at 
asymptotically large values of the chemical potential l30l . There, the mass of the diquark was defined as the energy of a bound state of 
two virtual quarks in the center of mass frame, i. e. as in our approach, in the rest frame for the whole diquark. 
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FIG. 4: The masses of diquarks. At /j. < /ic = 342 MeV, six diquark states are splitted into a (color)triplet of heavy states 
with mass Ma* and an anti-triplet of light states with mass Ma (see II55I I). In the 2SC phase (/i > Hc), one observes the 
Numbu-Goldstone boson, four light diquarks with the mass 3|/Lt8|, and a heavy singlet state with the mass M (solid line). The 
shaded rectangular displays the width F of the heavy singlet resonance; its upper border is half-width higher than the mass 
and the bottom border is half-width lower. 



where 

r ( 2 s f _A 1 f _fq 1 

"^^"^ J {2,rfEXmlY-pl]^ J (27r)3iJ-[4(i?-)2_p2]' ^^^^ 
J (2. [ d\ E+ r d^q E- 

J i2n)^E+[4{EX)^-pl] J i2n)^E^[4{E^r-piy ^ ' 

The mass spectrum is defined by the equation 

detr2«(po) s ipliipl - 4A2)/2(p2) „ 4/2(^2)1 ^ o_ (52) 

In the pQ-planc tliis equation has an evident zero, corresponding to a Nambu-Goldstonc boson, = 0. Detailed 
investigation, similar to that from 0, shows that on the second Riemann sheet of Pq there is another zero of H52|) 
which corresponds to a heavy resonance. Its mass M and width F are depicted in the Fig. 4. 

As a result, we conclude that in the 2SC phase there are four light real scalar diquark excitations with mass 
3|/xg|, and one Nambu-Goldstone boson, which appears due to a spontaneous breaking of the SUc(2) x 11^8(1) color 
symmetry down to SUc(2). Moreover, a heavy scalar diquark resonance, which is an SUc(2)-singlet, is also presented 
in the mass spectrum of the model at /i > /ic = 342 MeV (see Fig. 4). 



B. Scalar diquarks in the normal phase (A — 0, fis ~ 0) 

In the SUc(3)-symmetric phase (fi < /ic) the diquark gap A is zero and the three complex diquark fields A^(x) 
{A — 2, 5, 7) are not mixed with other fields in the second order effective action of the model. Moreover, at A = 0, as 
it is easily seen from the color neutrality constraint H22|) . we have /xg = 0. So, in order to study the diquark masses, 
it is enough to consider, e.g., the A§-diquark sector. In this phase the determinant of the inverse propagator matrix 
F^*(po) looks like (we use the rest frame, where p = 0) 

detr'^po) = r5f^.A|bo)ri^5Ai-bo) - r%,^.{po)T%,^.{-po), (53) 

where r^s.^^(po) is presented in l|45|l . (Note, the last equality in formula (|48|l for detF^*(po) in the 2SC phase is 
based on the usage of a non trivial solution A 7^ of the gap equation (|20|l . So, in the normal SUc(3)-symmetric 
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phase, where A = 0, it is not valid.) Taking into account the relation m> fi that is realized in the normal phase only 
(see Fig. 1), we see that E = ^Jq^ -\-w? > ^ and, as a consequence, i?* > are fulfilled in this phase. So, one can 
easily integrate in H45|l over go and obtain the following expression that is suitable only for SUc(3)-synimetric phase: 



where e = (po + 2^)^. Clearly, the diquark mass spectrum is defined by the equation detr^''(po) = 0, or by zeros of 
(|54|l . where the function Fs{t) is analytical in the whole complex e-plane, except for the cut Arn? < e along the real 
axis. (In general, Fs{e) is defined on a complex Riemann surface which is to be described by several sheets. However, a 
direct numerical computation based on eq. (|54|l gives its values on the first sheet only (we use the parameter set (|23f) '). 
To find a value on the rest of the Riemann surface, a special procedure of continuation is needed.) The numerical 
analysis of (|54|l on the first Riemann sheet shows that the equation r^s.^^(po) = has a root (cq) on the real axis 
(0 < £0 < 4r7i2), providing us with the following massive diquark modes which are the solutions of the eq. (|53|l 

(MJ,)2 = (1.998m- 2^)2, {Ml,f^{l.998m + 2fif. (55) 

We relate in (|55ll to the mass of the diquark with the baryon number B — 2/3 and M^, to the mass of the 
antidiquark with B = —2/3. (Qualitatively, a similar behavior of diquark and antidiquark masses vs. ^ was obtained 
in [sJl in the NJL model with two-colored quarks.) The difference between diquark and antidiquark masses in H55|) is 
explained by the absence of a charge conjugation symmetry in the presence of a chemical potential. 

Finally, due to the underling color SU(3)c symmetry, the previous statement is valid also for Ag, A5 and Ay, A7. 
As a result, we have a color antitriplet of diquarks with the mass (|55|l as well as a color triplet of antidiquarks 
with the mass M^, . The results of numerical computations are presented in Fig. 4 for ^ < Hc = 342 MeV. 

Recall that in our analysis, we have used the constraint Hg — 3G/4, thereby fixing the constant Hg through G. It 
is useful, however, to discuss now the influence of Hg on diquark masses. Indeed, it is clear from H54|l that the root 
eo lies inside the interval < eo < Arn^ only if H* < Hg < H** , where H* and H** are defined by 



4F,(4m2) 4 [AVm2 + A^ + m2 ln((A + Vm^ + A2)/m)] ' 

H** = 1 ^ 7r2 ^ 3mG 

' ~ 4F,(0) 4 [AVm2 + A2 - m2 ln((A + Vm^ + A^)/m)] 2(m-mo)' ^ ^ 

In this case, there are stable diquarks and antidiquarks in the color symmetric phase. The behavior of their masses 
qualitatively resembles that given by eqs. (|55(l . For a rather weak interaction in the diquark channel [Hg < H*), 
eo runs onto the second Riemann sheet, and unstable diquark modes (resonances) appear. Unlike this, a sufhciently 
strong interaction in the diquark channel {Hg > H**) pushes eo towards the negative semi-axis, i. e. {pq + 2fi)^ < 0. 
The latter indicates a tachyon singularity in the diquark propagator, evidencing that the SUc(3)-color symmetric 
ground state is not stable. Indeed, at a very large Hg, as it has been shown in [33| at Hp = 0, the color symmetry is 
spontaneously broken even at a vanishing chemical potential. We guess that this result remains true at rather small 
values of Hp (Hp < Hg) as well, justifying the above mentioned tachyon singularity of the diquark propagator. 



VI. MASSES OF PSEUDO-SCALAR DIQUARKS 

It is clear from (|35f) that in the framework of the NJL model ^ pseudo-scalar diquarks are not mixing with each 
other as well as with meson- and scalar diquark fields. So, to get their masses, we will start from the general expression 
(I41|l for the inverse propagator matrices of pseudo-scalar diquarks. 



A. Pseudo-scalar diquarks in the 2SC phase (A 7^ 0, 7^ 0) 

In the 2SC phase there is an SU(2)c symmetry between A|-A|* and A7-A7* sectors, so it is enough to study an 
inverse propagator, e.g., in one of these sectors. For A — 5 it follows from 141|) that 

where P(x),Q{x) = Ag(x), Ag*(x). Further, using the evident expression H37|l for iS^^gg, one can obtain the matrix 
elements (|57ll in the rest frame, p = 0, of the momentum space representation 



(E- -po + E+)0{-E-) E+-E+ 



^t^Al^P^) AHp^-J (2^)3 \ E^ L (p, _ E-Y - {ElY 2{p„ - E- - EX)E 
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iE+ + po + E-)e{E+) 



E- +E-. 



{po + E+f- 
E+-Ei 



2ipo + E+^E+)E+ 



2{po + E+ - E^)E^ 

{E- -po-E-)e{-E- 

{po-E-y-{E^r 



<f UE+ - E+ ~ po)9iE+) 

En + ^+)2 _ (^+)2 

E- +E^ 



2{po - E- - EX)E^- 



(58) 



r^/p .p, (p) 



r 



Ar A?(~P)' ^^'^ ^%A?(P) - -'"a?* A? 



0. Let ^ 



JJJ- + V- Then, using for jjj- the gap equation 



^P^P'\t^j — - Al'Al^ ^'^ """" " APAP^f ~ Af*Af*(P) 

(|20|l (recaU, in the 2SC phase A 7^ 0), we obtain from H58|) (transforming the multipher before the second square 
bracket in (inHJ as g*^/^^ ^ 1 - m^/E'^): 



T%..p{po) =ri + 2{po - 3^s)H{po) + m^H{pQ), 



(59) 



where H(jjq) is defined in H47|) . Since in the 2SC phase m << A (or m even a zero if mo = 0), we ignore the last 
term in H59|l (due to this reason an exphcit form of H{po) is not presented here) and obtain in this way the following 
expression for the determinant of the matrix T'^P{po)'. 



detr^P(po) - T%„^p{po)r%^p,{po) 



-p5p 
^ A?' A 



?(po)r^^p-Af("Po) 



7? + 2(po - 3/^8)i7(po) V-2iPo + 3^^s)H{-po) « 7^ + 2{po - S^isf ^ r) + 2{po + 2^JL^fp 



(60) 



where (3 = dH{pQ)/dpo\pg^3fj_g. In the last relation in H6Q(I we have expanded the function H{pq) into a Taylor-series 
of pq at the point pq = 3/i8, and took into account that H{3fj,s) equals zero under the color neutrality constraint (see 
the end of section lV^ . Solving in this approximation the equation detr^P(po) = 0, it is possible to find pseudo-scalar 
diquark excitations with two different masses: 



A'f - 



3/i8 + v-v/m 



(61) 



Since /3 < (as is easily seen from (07J), we conclude from formulae H61|) that at ry > both and are real 
(positive) quantities, suggesting that these pseudo-scalars are stable particles at Hp < Hg. (The case 77 < 0, which 
corresponds to an unstable 2SC ground state, will be discussed in details below, after (j^SJl-) The similar is true in 
the Ay- Ay* sector of the model, so in the whole A|-A|*, Ay- Ay* sector of the NJL model which is under the color 
neutrality constraint there are four massive pseudoscalar excitations: two of them form an SUc(2)-antidoublet with 
mass M^^, another two particles form an SUc(2)-antidoublet with mass 



D2- 



Now let us consider the 2SC ground state excitations in the Aj-Aj* sector. In this case, starting from the effective 
action H39|) . it is possible to obtain the inverse propagator matrix T^p which is defined by the relation (|41|l . In the 
rest frame, where p = {po-, 0, 0, 0), its Fourier-transformed matrix elements look like 



A^A: 



■{Pa 



1 



AH, 



p 



d^q r 2m2 Po + 90 - E"' 



{2n) 



qo 



E- 



E^ 



Po + qo 



E+ 



E^ 

qo + E^ 



{Po 



qoY-{ElY 
Po + qo + E' 



1l 



{ElY 

qo-E- 



(po + qoY 



{EX? 



I'o 



(Eir 



r 



2p 



'o) = 



d^q 



2to2 



iPo+qor-{E^r 

Tl^.Jvo) = r^'' 



^0 



{E-^Y 



(62) 



A^Af 



E-^ 



i(27r^4 
1 



AS* a: 
1 



p,{po) = -4A^P(po), 



Y I E^ (po + qo? - {EX? 



qo) 
1 



ql-iE-^Y 
1 



{po + 90)^ 



{eX)^ 



Ql 



{E^? 



(63) 



{E+r iPo + qo? - {E^ 

Using in (|62(l - (|63|l the substitution jjj- = jjj- + -q and then eliminating the coupling constant Hs in favor of another 
model parameters (with the help of the gap equation (|2U|) '). wc obtain after qo'integrations: 



P(Po) =/o(Po)+™'^(Po), 



n2p 
- ArA 



p{Po) (4A2 - 2pI)I^{pI) + Apohipl) + m''B{po), 



(64) 



where /o(po) and Ii{pl) are presented in (|50|l and (|51|) . respectively. The last terms in each of expressions (|64f) 
are proportional to m^. In the 2SC phase the constituent quark mass m is a vanishingly small quantity (or it is 
exactly zero if the current quark mass vanishes, mo = 0) as compared with A, etc (sec Fig. 1), so we will ignore the 
contributions of these terms in the matrix elements (|62|I - H63|I . Thus, there is no need to have explicit expressions for 
the functions A{pq) and B{pq) from H64|l . In this approximation the determinant of the inverse propagator matrix 
r2p(po) looks hke: 



detr2p(po) = r'^" 



^p.Ap(po)r^^p^p.(po) - rrPA;(^^o)r^p.^p 



n2p 



n2p 



{Po) 



i^-2pih{pi) ii-2{pi-A^')h{pi) -IGpil'M) 



(65) 
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Obviously, at 77 = the equation detT'^P{po) = coincides with eq. (|52|l and has the same solutions. The first one, 
Pq = 0, corresponds to a stable massless pseudo-scalar excitation, the second one lies in the second Riemann sheet of 
the variable Pq. So it is a heavy pseudo-scalar resonance and its mass and width are represented in the Fig. 4. At 
small nonzero values of 77 it is reasonable to suppose that the equation detT'^P{po) — has a root lying on the second 
Riemann sheet of pg as well. It might be considered as a weak disturbance of the resonance solution of this equation 
at 7y = 0, so its mass and width behavior vs ^ are qualitatively the same as on the Fig. 4. Another solution of this 
equation, Pq = (Afj^g)^, should not be significantly different from the solution Pq = at r/ = 0. So, in searching of 
(Mj^g)^, one can expand the expression into a Taylor-series of p§ = 0: 



1662 -I- 2a{Tj + 8A2a) + 2r]{a - iA^a' 



(MP )2 '/V'/ ^ °^ (66) 



where a = /o(0), b = /i(0), a' = Iq{0). Note that both the heavy resonance and stable excitation with mass squared 
in the pseudo-scalar diquark channel are singlets with respect to SUc(2). Since a > 0, the expression is a 
positive one at rather small and positive values of ry. However, at sufficiently small, but negative values of rj, it is a 
negative quantity, i.e. a tachyonic pseudo-diquark excitation appeared in the model. This fact indicates the instability 
of the 2SC ground state with (A|(a;)) = A, (A§(a;)) = 0, (Af (a;)) = 0, and (A^(a;)) = 0. Perhaps, in this case, i.e. at 
Hp > Hs, the phase with nonzero ground state expectation values of pseudo-scalar diquarks, (A^(x)) 7^ 0, should be 
realized. 

As a result, we have shown that in the 2SC phase of the NJL model there are five stable diquark excitations 
in the pseudo-scalar channel. They form a singlet as well as two antidoublets (in the case Hp < Hg) of the SUc(2) 
group with masses presented in H66|l and H61|l . correspondingly. Moreover, there is also a heavy resonance that is 
SUc(2)-singlet with mass about 1100 MeV in this channel. 



B. The case of normal phase (A — 0, fig — 0, m 0) 

Suppose that we are in the SUc(3)-symmetric (normal) phase of our model, where A = 0. As it is easily seen from 
the color neutrality constraint H22() . in this phase fig = 0. Moreover, here = = E ± ij, > 0, since in this phase 
m > fj, (see Fig. 1 at < /i^). Then, without loss of generality, it is sufficient to study the mass spectrum, e.g., in 
the sector of Ag-Ag* diquarks. For the normal phase we have from H58|) : 

- - ^^J^'E4E^-{l + 2^y ^m-p- ^^^'^^ 

where e = (po + 2/^)2 and the function Fp{e) is increasing on the interval {—00, Am?). Moreover, it is analytical in 
the whole complex e-plane, except for the cut 4m2 < e along the real axis. The masses of pseudo-scalar diquarks are 
defined by the equation 

detr5^'(po) = r^f/.^, {po)T%^., ipo) - r%,^, iPo)T%,^. i-po) - 0, (68) 

5 5 5 5 5 5 5 5 

i.e. by zeros eo of the matrix element (|57|l such that < eq < 4m2 or lying in the second Riemann sheet of the complex 
variable e. (The first one correspond to masses of stable excitations, the second one - to masses of resonances.) In 
the present consideration we restrict ourselves to looking only for stable pseudo-scalar diquarks. It is clear that there 
is a single zero eo of (|67|l . obeying the condition < eo < 4to2^ if and only if the coupling constant Hp is constrained 
by the relation H* < Hp < H** , where 

h: ' 



P 4Fp(4m2) 4[A%^?TA2-m2ln((A + %A^?TA2)/m)] 



1 7r2A\/"i2 + A2 , 

P " 4Fp(0) " 4 [3to2A2 + A4 - 3m2AVm2 + A2 ln((A + Vm^ + A^)/m)] ' 



(Note, Hp — H** from H56|l.) In this case the masses of stable pseudo-scalar diquarks and antidiquarks are the 
following 

(Mi)' = (^/^-2M)^ (MX.)'-(V^+2/.)^ (70) 

respectively. (The mass splitting in (|70|l is again explained by the absence of a charge conjugation symmetry in the 
presence of a chemical potential.) It follows from the underling color SU(3)c symmetry of the normal phase that at 
Hp < Hp < Hp* there is indeed a color antitriplet of pseudo-scalar diquarks with the mass M£ as well as a color triplet 
of pseudo-scalar antidiquarks with the mass M^, (see (|70|l '). For other regions of the iJp-values, stable pseudo-scalar 
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diquark excitations of the SUc(3)-color symmetric ground state are forbidden. For a rather weak interaction in this 
channel {Hp < H*), eo runs onto the second Riemann sheet, and unstable pseudo-scalar diquark modes (resonances) 
appear. Unlike this, a sufBciently strong interaction in this channel {Hp > H**) pushes eo towards the negative 
semi-axis, i. e. {po + 2fj,)'^ < 0. The later indicates a tachyonic singularity in the pseudo-scalar diquark propagator, 
evidencing that the SU(3)c color symmetric ground state is not stable. In this case we guess that a parity-breaking 
color superconducting phase is realized in which the ground state expectation values of pseudo-scalar diquarks are 
not zero, (A^(x)) ^ 0. 

Finally, a few words about diquarks in the particular case (recall, /i < fic), Hs — Hp = H, which corresponds to a 
NJL model inspired by a one-gluon exchange approximation in QCD. In this case we see that ii < H < H*, then 
both scalar- and pseudo-scalar diquarks are resonances. If H* < H < H**, then in the normal phase, including the 
case fj, = 0, only scalar diquarks are stable, but pseudo-scalar ones are unstable particles. For larger values of H the 
normal phase is unstable in itself, since either or both scalar diquark propagator (at H** = H* < H < H**) and 
pseudo-scalar one (at H** < H) have tachyonic singularities. Hence, one can conclude that at Hg — Hp = H scalar 
diquarks are allowed to exist, at H* < H < H** , as stable excitations of the normal phase. Pseudo-scalar diquarks 
in this phase are always unstable particles (resonances). 

VII. SUMMARY AND DISCUSSION 

In our previous papers 0,0) the masses of mesons and diquarks, surrounded by moderately dense quark matter, 
were investigated in the framework of NJL model ^ at Hp = 0, and the color neutrality constraint was missed, 
for simplicity. In the present paper, we have calculated the mass spectrum of meson and diquark excitations in the 
color neutral cold dense quark matter. We started from a low-energy Nambu-Jona-Lasinio type effective model 
for quarks of two flavors, with a quark chemical potential fi and extended by including the chemical potential /is 
of the 8th color charge. Moreover, the interaction in the pseudo-scalar diquark channel was taken into account, in 
addition. We considered only the interplay between normal- and 2SC phases. This is a quite reasonable assumption 
in the framework of the model . Then, it was shown that in the presence of color neutrality the transition to the 
2SC phase occurred at a somewhat smaller value of the quark chemical potential (/ic — 342 MeV) than without this 
constraint (/ic = 350 MeV). 

It was proved in the present paper that in both models and i.e. with or without a color neutrality constraint, 
the (T-meson is mixed with the scalar diquark A| in the 2SC phase. In the previous paper this mixing was ignored 
in the consideration of the cr-meson mass. At first, we have found that, if cr-A| mixing is ignored as in 9], then 
the color neutrality requirement does not change qualitatively the properties of tt- and cr-mesons, obtained in the 
framework of NJL model without the /ig term. This is an expected result, since both models and ^ have 
an identical chiral symmetry. Hence, at small values of /ig (see Fig. 2) the meson masses aquire small corrections as 
well (compare Fig. 2 from ^ and Fig. 3 of our present paper). It follows from our consideration that in this case 
(without mixing) both a- and tt- meson are stable particles in the 2SC phase with masses of about 340 Mev (Fig. 3). 
However, if mixing is taken into account, then in the 2SC phase the u-meson is a resonance, decaying into a pair of 
quarks with a rather small width ~ 30 MeV (see Appendix). As far as we know, the properties of tt- and cr-mesons 
in the 2SC phase have not been discussed in the literature before. 

Moreover, the properties of scalar diquarks in the 2SC phase are changed drastically, when the color neutrality 
condition is imposed. Indeed, for the model (Q) we have found in the 2SC phase an anomalous number of three 
Nambu-Goldstone (NG) bosons, the SUc(2)-antidoublet of light diquarks, and a heavy resonance that is an SUc(2)- 
singlet H,I3- Contrary, our present investigation shows that in the model Q the scalar diquark sector of the 2SC 
phase contains two real SUc(2)-antidoublets of light excitations with the same mass 3|/i8|, one Nambu-Goldstone 
boson, and a heavy resonance with mass about 1100 MeV (see Fig. 4). To understand such a sharp difference in 
scalar diquark masses, predicted by these two models, it is necessary to compare their color symmetries. The first 
model, Lagrangian is invariant under SUc(3). However, in the second model (01 this symmetry is broken explicitly, 
due to the presence of the /zg-term, to the subgroup SUc(2) x Ua8(1). Then, in the 2SC phase, where the ground 
state is an SUc(2) invariant for both models (here (A|(x)) ^ 0, (A| j{x)) = 0), we have spontaneous breaking of the 
above mentioned symmetries. As a consequence, there are five broken symmetry generators and an abnormal number 
of three NG bosons for the model ^ (the explanation of this fact is presented in details in [81). On the other hand, 
for the model l(2Jl we have in the 2SC phase only one broken C/^g-symmetry generator, resulting in a single NG boson. 

The properties of the pseudo-scalar diquarks in the 2SC phase depend essentially on the relation between coupling 
constants Hs and Hp. First, note that at Hp > Hs the 2SC phase is an unstable one due to a negative mass squared 
(|66|l of an SUc(2)-singlet pseudo-scalar mode (tachyonic instability). At Hp < Hg the pseudo-scalar excitations of 
this channel form in the 2SC phase two real stable SUc(2)-antidoublets with different masses (|til|) as well as the stable 
light SUc(2)-singlet with mass and a heavy resonance with mass ~ 1100 MeV. 

We have also found that the antidiquark masses exceed those of the diquarks in the normal SUc(3) symmetric phase 
(for n < 342 MeV). This splitting of the masses is explained by the violation of C-parity (charge conjugation) in 
the presence of a chemical potential. In contrast, at /it = the model is C-invariant and all diquarks and antidiquarks 
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of the same parity have an identical mass. It follows from our investigation that stable quark pair formation occurs 
in the scalar channel at a weaker coupling strength than in the pseudo-scalar one. Indeed, if the parameter set of the 
model is fixed by the relations 123|l , then in the normal phase scalar diquarks are stable particles (with masses about 
700 MeV at ^ = 0), since in this case H* < Hs = 30/4 < H** (see (f^ ). However, only at sufficiently high values 
of Hp, i.e. at H** = H* < Hp < H** (see (I69ll ). stable pseudo-scalar diquarks might exist in the normal phase. On 
the other hand, if Hp < Hs = 3G/4, then in the normal phase pseudo-scalar diquarks are resonances, decaying into 
two quarks. If ^ exceeds the critical value and the system passes to the 2SC phase at Hp < Hs = 3G/4, then five 
of six pseudo-scalar diquarks aquire stability. Really, in the 2SC phase, in contrast to the normal one, all particles 
move inside the medium, so their decay might be prohibited by a Pauli blocking principle (Mott effect). Suppose 
that Hp = Hg = H (in this case the NJL model is inspired by a one-gluon exchange approximation in QCD). Then 
at H < H** the normal phase is a stable one, but at H > H** a tachyonic instability appears, so the normal phase 
is destroyed (see the section |VB|) . We have shown for this particular case that in the normal phase there might exist 
stable scalar diquarks, however pseudo-scalar diquark modes are always unstable excitations. 

Of course, all observable particles render themselves as colorless objects in the hadronic phase, and the diquarks 
are expected to be confined, as they are no SU(3)c color singlets. Nevertheless, one may look at our and other 
related results on diquark masses as an indication of the existence of rather strong quark-quark correlations inside 
baryons, which mig ht help to explain baryon dynamics. Some lattice simulations reveal strong attraction in the 
diquark channel |33j| with a diquark mass ~600 MeV. Recently, in |34j , the mass and extremely narrow width, as well 
as other properties, of the pentaquark 0+ were explained just on the assumption that it is composed of an antiquark 
and two highly correlated ud-pairs. At present time, the nature of the mechanism which may entail strong attraction 
of quarks in diquark channels is actively discussed both in the nonperturbative QCD and in other models (see, e.g., 
[35I I and references therein). 

For simplicity, we have studied the masses of the one-particle excitations of the color neutral quark matter. Similar 
investigations can be performed both in the color neutral and /3-equilibrated two-flavor NJL model, where electric 
charge neutrality is required in addition, so one more (electric charge) chemical potential /ig must be introduced. In 
this case, for some range of the coupling constant Hs the gapless 2SC phase (g2SC) is realized (see, e.g., 3, where the 
particular case Hp = of the NJL model ^ was considered). In contrast to the ordinary 2SC phase, two additional 
gapless quark excitations then appeared in the g2SC phase, but its meson and diquark mass spectrum, presumably, 
will not changed qualitatively. Our belief is based on the structure of the ground state expectation values of scalar 
diquarks in the g2SC phase, i.e. (A|(a;)) = A 7^ 0, (A| 7(2;)) = 0. This is expected to be identical to that of our 
present consideration when only the color neutrality of the quark matter is taken into account. 
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APPENDIX A: CORRECTION TO THE a-MESON MASS 



In section HVI the numerical values for the ct- meson mass was obtained in the assumption that the mixing of a 
with A| is absent. Moreover, we have shown in this approach that the tr-meson is a stable particle in the 2SC phase. 
Now let us prove that a deeper investigation of the tr-meson mass, based on the inclusion of the mixing between a 
and A2, results in the conclusion that in the 2SC phase cr is a resonance, having a small decay width into a quark 
pair. 

Indeed, starting from the total number of effective actions H18f) . (|28ll . and H38|) it is possible to find the 3x3 inverse 
propagator matrix Q^^ of a-, A|-, and A|* fields. In the center of mass frame of the momentum space representation 
{p = 0) it has the following form in the 2SC phase: 

r(po) , r<^A|(po), r<^A|-(po) 
g-'ipo)^ { ^AiApo), rj^s^^^ipo), ri^.^„(po) |, (ai) 
TArApo), riVAsbo), riVArbo) 



where 



"<tA|(Po) = TA-*a{Pa) = ^aA-*{-Po) = ^A^ai-Po) = 

= 4 A f [ ^E++po 2E--P0 

J (27r)3 \EE+[AiE+r ~ p^] ^ EE^[A{E^r - p^,] ^ ' ^ ^ 
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Other matrix elements from (jAlj) are presented in (|32|l and (|49|l . Note that the matrix elements ljA2p . mixing cr with 
A2 and A2*, are proportional to the dynamical quark mass m. So, in the 2SC phase both m and these matrix elements 
may be considered as small quantities (see Fig. 1). The mass spectrum is defined by the equation det{G~^{po)) = 0, 
which has a rather complicated form (note, det(^~^(po)) is an even fimction vs po, i.e. it depends on p^): 

FipDToipD^m^A'fipl), (A3) 

where 

- r,A.rA-.ri''^Aj - r.A*rA,.ri''.^^ - f{pI)t,{pI), (A4) 

and the functions Fq, Fi are defined in H33|) . (|34l) . In the pp-plane the equation ljA3|l has three solutions. One of them 
corresponds to a Nambu-Goldstone boson, pg = 0. The other two we denote as {pq)a and (Po)cr- Since m is a small 
quantity, the right hand side (RHS) of (|A3I) can be considered as a small perturbation. Then the solution (p§)cr (as 
well as (po)a) can be constructed in the framework of a perturbative expansion, based on the smallness of the RHS 
of ljA3|l . In the zeroth order we have F{p^)rQ{pl) = 0, so = M^, where M„ is the zero of the function Fo(pq) 

and is graphically represented in Fig. 3. One can easily obtain the first perturbative correction 

(^°^-^^'- + f(Af^)Fg(A4^) +--- 

and so on. Now let us put our attention on the fact that /(pg) and F{pq) are analytical functions in the whole 
complex pg-plane except the cut, composed from all real points such that 4A^ < Po < oo- (In the rest of the real 
Pg-axis /(po) and F{pq) take real values.) Through the cut these functions can be analytically continued to the second 
Riemann sheet. In our case Mo- « 350 MeV and A ^ 100 MeV, so G (4A^,oo), i.e. it lies on the cut. Hence, 
both f{M^) and f{M^) have imaginary parts, and the solution IjASp lies in the second Riemann sheet. It means that 
(pg)a- corresponds to a resonance in the mass spectrum. Numerical estimates show that the width of this resonance 
is a rather small quantity, less than 30 MeV. 

Similar corrections can be easily performed for the heavy scalar diquark resonance mass, corresponding to another 
solution (po)a of the equation HA3p . 
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